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Abstract 


The dynamic behavior of a five cells proton exchange membrane fuel cell (PEMFC) stack operating in dead-end mode has been studied at room 
temperature, both experimentally and by simulation. Its performances in “fresh” and “aged” state have been compared. The cells exhibited two 
different response times: the first one at about 40 ms, corresponding to the time needed to charge the double-layer capacitance, and the second one 
at about 15—20s. The first time response was not affected by the ageing process, despite the decrease of the performances, while the second one 
was. Our simulations indicated that a high amount of liquid water was present in the stack, even in “fresh” state. This liquid water is at the origin 
of the performances decrease with ageing, due to its effect on decreasing the actual GDL porosity that in turn cause the starving of the active layer 
with oxygen. As a consequence, it appears that water management issue in a fuel cell operating in dead-end mode at room temperature mainly 
consists in avoiding pore flooding instead of providing enough water to maintain membrane conductivity. 
© 2007 Elsevier B.V. All rights reserved. 


Keywords: Proton exchange membrane fuel cell; Water management; Dead-end mode; Stack performance; Startup; Dynamic behavior 


1. Introduction At the anode side, hydrogen gas which acts as fuel is oxidized 
according to the following reaction: 

As schematically fuel cell assembly shown in Fig. 1, a proton 
exchange membrane fuel cell (PEMFC) consists of two elec- 
trodes containing a thin layer of catalyst in contact with a plastic 
membrane separating gas supply chambers. The electrodes are 
constituted by a gas diffusion layer (GDL) and an active layer 
(AL). The GDL is usually made of carbon cloth, carbon powder O2 +4H* +4e7 > 2H20 (2) 
and PTFE, while AL usually consists in a Pt/C powder uni- 
formly agglomerated with Nafion® ionomer. Both layers have 
a porous structure. This assembly is sandwiched between two 
current collectors [1]. 


Hz —> 2H* +2e7 (1) 


Similarly, oxygen enters the fuel cell at the cathode where it 
is reduced in water according to the following reaction: 


Gaseous oxygen fed in the gas channel diffuses through the 
GDL and in the pores of the AL, where it dissolves in Nafion® 
on the surface of the agglomerates. Dissolved oxygen diffuses 
towards the centre and reacts on the numerous Pt particles. The 
management of the water generated by reaction (2) remains a 
~~ major issue. On one hand, its shortage causes a dehydration of 
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Nomenclature 

{] concentration (mol m~?) 

C capacitance (F m7?) 

D diffusion coefficient (m? s7!) 

F Faraday constant (C mol`!) 

j current density (A m7?) 

M molar mass (g mol`!) 

n number of moles 

P pressure (Pa) 

R resistance (Q) 

R Gas constant (J K7! mol!) 

sS volume fraction of water in the GDL 

S geometrical surface area (m?) 

t time (s) 

T temperature (K) 

U cell voltage (V) 

V volume (m?) 

x distance in direction orthogonal to electrode sur- 
faces (m) 

Greek symbols 

ô GDL thickness (m) 

E porosity 

n overvoltage (V) 

p volumetric mass (kg m7?) 

Subscripts 

0 initial 

act reactive 

AL active layer 

C capacitive 

dl double layer 

dry in dry state 

e electrode 

eff effective 

eq equilibrium 

exp experimental 

F Faradaic 

GDL gas diffusion layer 

int internal 

imp imposed 

sat saturation 


for the improvement of the water management, since it allows a 
better understanding of its effect on fuel cell performances. 

The cell would be operated on 100% hydrogen and oxygen 
utilization, i.e. on a stoichiometry of 1. This is also referred 
to as “dead-end” operating mode. From the system point of 
view, dead-end operation can be accomplished by closing the 
gas outlet of the cell. 

In the past, the mathematical models of the PEMFC were 
usually restricted to the steady state conditions and transient 
models only recently appeared. But these models are so far 
restricted to stoichiometric mode. Bernardi [2] first proposed 
a one-dimensional model in order to study the water man- 
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Fig. 1. Schematic representation of a PEMFC and the gas diffusion electrode. 


agement and to identify the humidification conditions which 
induce either the dehydration of the membrane or excessive 
flooding. The sensitivity to the water balance of the PEMFC per- 
formance was demonstrated. Some other models were derived 
from the original one to take into account heat management 
[3], mass transport in the gas diffusion electrode [4] or to 
introduce a different treatment of the electrochemical reaction 
[5]. Bernardi et al. [6] and Springer et al. [7] also presented 
a one-dimensional model to investigate the factors that limit 
cell performance and to elucidate the mass transport mecha- 
nism within the complex network of gas, liquid and solid phases 
constituting the gas diffusion electrode. Contrary to the previ- 
ous models, this latter one considered the molar changes along 
the gas channel. These approaches made it possible to eval- 
uate the losses in the cathode gas backing and the catalyst 
layer. 

So far, efforts have been made to model the problem 
in three dimensions. A higher dimensionality makes it pos- 
sible to describe the hydrodynamics and multi-component 
transport inside the flow channel for reactant distribution. Multi- 
component models were thus developed for cathode or for whole 
PEMEC, for both single-phase [8] and two phase flow [9]. These 
models showed that the gas distributor geometry has a signif- 
icant effect on the diffusion of the reactants and products in 
determining the performance of the cell. Moreover, they con- 
firmed that the performance of cathode is strongly influenced by 
the presence of liquid water and its removal rate, especially at 
high current density. 

The dynamic study of a fuel cell appears to be of great interest 
to provide detailed understanding of mass and charge transport 
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through the gas diffusion electrodes and is of extreme impor- 
tance for the control strategy and system management in power 
generation systems, especially when there are power injections 
into the network. Thus, a first equivalent circuit was proposed to 
simulate the impedance spectrum of a PEMFC [10]. This study 
showed that the impedance of a fuel cell is a powerful tool in 
order to characterize the intra-electrode processes occurring in 
gas diffusion electrodes. 

Many recent dynamic models of PEMFC [9,1 1-19] have been 
developed, based on the physical and chemical knowledge of 
the phenomena occurring inside the cell. Some of them even 
take into account the multi-dimensional geometry of the cell 
[9,18,19] while some others exhibited an important effect of 
water management [13-15]. Ceraolo et al. [11] provided the 
static and dynamic model of PEMFC. Their one-dimensional 
model only considered the cathode side, and the transient terms 
were added to the species balance and phase potential equa- 
tion. They obtained a good correlation between their numerical 
predictions and experimental results. Yerremalla et al. [12] pro- 
posed a dynamic model based on semi-empirical equations. 
In contrast, the model of Wang and Wang [14] includes the 
full description of water and proton co-transport in the three- 
dimensional MEA, thus it could provide detailed information 
of water behaviors in this vital component, such as water accu- 
mulation. It is found that membrane hydration occurs over a 
period of 10s, the gas transport of 0.01—0.1 s, with the double- 
layer discharging being negligibly fast. Some other models were 
developed in a more “System” approach [20]. These models are 
relatively simple which allow correctly simulating the behav- 
ior of the fuel cell inserted in an electrical network. The main 
disadvantage of these models is that they are far away from the 
physics of the fuel cell, i.e. their parameters are physically non- 
representative. In parallel, some other multi-scale models were 
elaborated [21] to predict the dynamic and steady state behav- 
iors for the triple contact (where the electrochemical reaction 
proceeds). These models made it possible to describe quanti- 
tatively the reaction mechanisms, the polarization curves and 
impedance spectra of the fuel cell. However, all these mod- 
els deal with a cell operating in stoichiometric mode and little 
attention has been paid to systems that operate in dead-end 
mode. 

In addition, few experimental studies of dynamic behavior of 
the PEMFC were performed, all related with operation in stoi- 
chiometric mode with changes between two operating regimes 
{11,13,22—24]. None of these works deals with the starting 
behavior of a PEMFC operating in dead-end mode or studied 
the transient behavior to a step current/voltage from open cir- 
cuit voltage (OCV). Kim et al. thus studied, using a two channel 
digital oscilloscope, the influence of stoichiometry [22], reser- 
voirs and fuel dilution [23], flow field design and voltage step 
magnitude [24] on the dynamic behavior of PEMFC operating 
in stoichiometric mode. Their results indicate that the response 
time to a step voltage was ranging between 3s for a change 
in operation from excess to normal stoichiometry to 15s for a 
change in operation from starved to normal stoichiometry [22]. 
All these studies, depending on the operating conditions, show 
either undershoot (transient current below the steady state value) 


or overshoot (transient current above the steady state value) 
behavior usually attributed to local changes in the polymer elec- 
trolyte conductivity caused by the change in local water content 
in the vicinity of the reactive sites [11]. 

In this paper, the authors have focused on the dynamic mod- 
eling of the transient response of a PEMFC stack operating in 
dead-end mode. So as to better understand this dynamic behav- 
ior, we studied both experimentally and numerically the starting 
performances of such cells. The experimental results were com- 
pared, in the first ten seconds of operation, with the prediction 
of a simple mono-dimensional model of the transient behavior. 
In addition, the model’s ability to predict the voltage level and 
the time response of the fuel cell was evaluated. 


2. PEMFC dynamic model 
2.1. Basic assumptions 


A dynamic one-dimensional PEMFC model is developed to 
predict the transient response of the PEMFC to a step current. So, 
a dynamic model of PEMFC capable of characterizing transient 
phenomena is proposed which incorporates simultaneously the 
prominent dynamic aspects occurring within the GDL and elec- 
trochemical double-layer charge or discharge. Based on the basic 
consideration reported in the literature, the main assumptions of 
the proposed model are as follows: 


(i) The influence of the anode side of the PEMFC on its tran- 
sient behavior was neglected, due to on one hand the faster 
kinetics of the hydrogen oxidation reaction with regard to 
the kinetics of oxygen reduction reaction [18], and on the 
other hand, to the faster diffusion of hydrogen gas in the 
anodic GDL with regard to the diffusion of oxygen on the 
cathodic side [25,26]. 

(ii) Due to the small time periods considered in this model, the 
system was considered isothermal. The study was focused 
on PEMFC starting performances at room temperature. 

(iii) The water vapor pressure in the gas phase was consid- 
ered as constant and fixed to its saturation value. This 
later hypothesis is reasonable as we consider a fuel cell 
operating at room temperature. The water produced by 
the electrochemical reaction was thus considered in liq- 
uid state and its effect was taken into account through an 
homogeneous decrease of the GDL porosity ¢, according 
to the following equation [27]: 


E = Edry(1 — 8) (3) 


where éqry and s, respectively, correspond to the GDL 
porosity in dry state and to the GDL water content (volume 
fraction of water inside the pores of the GDL). 

In addition, we also assumed that this liquid water partly 
blocks the electrode geometrical surface area Se, reducing 
the reactive area Sact with a magnitude proportional to its 
volume fraction, according to the following relation: 


Sact = Se(1 — s) (4) 
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(iv) The gas transport in the GDL was assumed to obey Fick’s 
law and an ideal gas mixture is considered. 

(v) Electrodes were considered as a planar interface at which 
reaction proceeds. The complex structure of the AL and 
the distribution of the reaction in its volume were not taken 
into account. 

(vi) The PEMFC was considered as operating in dead-end 
mode, at a fixed gas pressure, with periodic short bleeds. 
But, due to the short period of time considered in this study, 
no bleed was taken into account in our model. 

(vii) Due to the small time periods considered in this model, the 
membrane water content was assumed to be constant and 
no electro-osmotic flux or water diffusion flux was taken 
into account. 


2.2. Gas mass balance equation 


The oxygen entering the cell before reaching the catalyst layer 
diffuses, through the gas diffusion layer, within a gaseous mix 
constituted by water vapor and oxygen itself. This phenomenon 
is described by the equation of continuity for oxygen gas: 


d[O2] d d[O] 
dt — dx (Pero: 121) =0 (5) 


where Do, eff is the effective oxygen diffusion coefficient in 
the porous medium, that can be determined from the oxygen 
diffusion coefficient in free medium Do, and from medium’s 
porosity £, according to the Bruggeman equation [28]: 


Do;,ett = Do,é'” (6) 
And according to ideal Gas law: 


NO» 
Po, = y Rt = [O2|RT (7) 


Therefore, Eq. (5) can be rewritten as follow: 


dPo, d Po, 
D 2 = 8 
dt eff, O2 dx2 ( ) 


This mass balance equation has to be solved in the respect of 
the following boundary conditions: 


(i) At gas channel/GDL interface, the oxygen pressure is 
constant equal to the cathode pressure minus the water 
saturation pressure: 


Po, (x = 0) = Po — Px50,sat (9) 


(ii) At electrode/GDL interface, the oxygen flow rate is directly 
related to the Faradaic current jp, according to the Faraday’s 
law: 


dPo RT 
—D — (x = 6) = ——— jF. 10 
O2,eff Gy (x = ô) ar iF (10) 


The Faradaic current can be related to the overvoltage at 
the interface, n, through an approximated formulation of the 


Butler—Volmer equation, already used in the literature [11,16]: 


, Po, (x = ô) (F/2RT)y 
= 2 -1 
=m (a =5j=0" 
_ [Po = 8) 
jea (G a aes 1) (11) 
OQU = ) 


In Eq. (11), the number of electron involved in the reaction is 
equal to 1, while value of 4 would have been expected according 
to reaction (2). However, Parthasarathy et al. [29] have shown 
that, whatever the temperature in the range 30-80 °C, the oxygen 
reduction reaction at platinum/Nafion® interface proceeds with 
two different mechanisms, depending on the applied current den- 
sity. For fuel cell current densities above 107? A cm, areaction 
rate in first order with respect to O2 pressure was observed 
indicating a single electron transfer reaction. This behavior was 
reported to be linked with a reaction that proceeds on a Pt-oxide 
free surface [29]. 


2.3. Charge balance equation at the cathode/GDL interface 


An electrostatic electrical equilibrium is established at a 
metal electrode/electrolyte (Nafion®) interface resulting in a so- 
called “double layer” of separated negative (electrons inside the 
metal catalyst) and positive charges (protons in the Nafion®) 
across the interface [30]. The interface behaves like a capac- 
itor in which the charge density of electrons and ions at the 
interface depends on the potential difference across this double 
layer. As a result, when a current is set to a cell, the operating 
voltage shows an immediate change due to the internal resis- 
tance, usually called “ohmic drop”, and then evolves slowly as 
double-layer capacitance (Cg) is charged/discharged and elec- 
trochemical reaction proceeds simultaneously. The set current 
is thus split in a capacitive current jc “used” to charge/discharge 
the double-layer capacitance and in a Faradaic current jp, “used” 
to perform the electrochemical reaction [30]. 

Moreover, the capacitive current is unaffected by the liquid 
water present in the gaseous pore since charges accumulate of 
both sides of platinum particles/ionomer interface which are not 
impacted by liquid water. Oppositely, the Faradaic current is 
affected by the flooding of the gas pores because the access of 
oxygen gas to the active particles is hindered. As a consequence, 
the charge balance equation can be written as follow: 


; . ; ; dn 
JimpSe = jC Se + jr Sact > Jimp = Cal a + jr — s) (12) 


where jimp is the imposed current density. 
The overvoltage evolution resulting from a step current can 
thus be calculated as follows: 


dn jimp — (l — s) (13) 
dt ~ Cal 


The operating voltage evolution following a current step can 
in turn be calculated as follow: 


U(t) = Vocy — Rint jimpse — nt) (14) 
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2.4. Water balance equation 


If all produced water is considered to condensate in the GDL, 
homogeneously reducing its porosity, the change of the water 
content in the GDL is expressed from the Faradaic current: 


ds ip S; 
say VopL > = 1 e (15) 


PH20 
My30 


2 


2.5. Initial conditions 


This approach has been developed to simulate the dynamic 
response of the cell potential to a steep change of the current 
density from 0 to 0.5 Acm™?. So, we focused on the starting 
behavior of the fuel cell. As a consequence, the initial condi- 
tions are zero current and voltage equals the OCV. When no 
current passes in the cell, the oxygen pressure is homogeneous 
and equals the pressure at x = 0, as defined by Eq. (9). The initial 
value of the liquid volume fraction, so, is one of the model’s 
parameter. 


2.6. Resolution 


The GDL was divided in 50 elements in x-direction along its 
thickness. Eq. (8) was then discretized using a finite difference 
method. The overvoltage and the water volume fraction changes 
are calculated according to Eqs. (13) and (15), respectively. 
There are three adjustable model parameters to fit experimental 
data: the initial water volume fraction sọ, the equilibrium current 
density jeq and the double-layer capacitance Ca. The model was 
solved in two successive steps with Matlab, using the ODE 15S 
algorithm with the set of constants and parameters reported in 
Table 1. 

Time constants for electrochemical double-layer, and gas 
transport were reported to occur over a period of 0.01-0.1s 
[14]. In the first step, we thus solved the set of Eqs. (3)-(14) 
with time scale of 0.1 s to determine the voltage evolution asso- 
ciated with oxygen mass transport and double-layer capacitance 
charge. During this short time period, water accumulation effect 
was neglected. The value of the liquid fraction was fixed at so. 


Table 1 

Model constants and parameters 

Constants Parameters 
Po =2 bars =2.027 x 10° Pa jeq (Am~*) 
Pt50,sat = 1704 Pa SO 
T=288K Ca (Fm?) 
Jimp = 5000 A m~? 

S=0.03 m? 


8=3.54 x 1074 m (given by the manufacturer) 
Eary = 0.8 (given by the manufacturer) 

Do, = 1.192 x 1075 m? s™! [25] 

Uocy = 1.05 V 

Rint,resh = 1.05 x 1077 Q 

Rintaged = 1.3 x 10-7 Q 

pmo = 1000gL7! 

Mpo = 18 gmol7! 


Inthe second step, we used the voltage and the pressure profile 
calculated in the first step as initial conditions. We solved the set 
of Eqs. (3)-(15), considering that the double-layer capacitance 
is fully charged and water accumulation was taken into account. 

The unknowns that the model determines are the pressure 
profile inside the GDL, the cell voltage and the water volume 
fraction. 


3. Experimental 


The dynamic behavior of a PEMFC stack consisting of five 
cells with 300 cm? active area, developed by Helion company 
was studied by setting a step current signal (107 A s~! rate) from 
0 to 150 A (0.5 Acm~?), at a temperature of 288 + 2 K. Analyt- 
ical grade hydrogen and oxygen were used in this study as fuel 
and oxidant. 

The fuel cell test bench used in this study was operated by 
a Schneider Electric Premium system that independently con- 
trolled fluidic and electric circuits and guarantee the safety of 
the experiment (Fig. 2). The stack was fed with dry reactive 
gases and operated in dead-end mode. The pressure in each 
compartment was automatically regulated at 2 bars. Current 
was controlled by an Agilent N3300A electronic load with two 
N3306A modules (0-120 A, 0-60 V, 600 W). Before applying 
the current, the stack was swept for at least 2 min with reac- 
tive gases and then maintained at operating pressure for a few 
minutes. When set, current was maintained during 30 s and then 
steeply shut down. The temperatures of inlet and outlet cells of 
the stack were measured by thermocouples during the whole 
experiment and did not increase more than 2 K. Current, total 
voltage and each cell voltage evolutions were recorded using 
a Keithley KUSB-3108 data acquisition module, at a 2000 Hz 
frequency in the first second after the current setup and then at 
a 5 Hz frequency from 1 to 30s after the current switch on. 

The stack dynamic behavior was first studied just after man- 
ufacturing and few hours of conditioning operations, and then 
after cycling and ageing during 3 months in conditions that can- 
not be detailed for confidentiality reasons. The former case will 


Hydrogen’ 
sensor 


Gases 
humidifyers 
Thermostatic bath: 
regulation 


K-USB Data 
acquisition 
module 


W Helion Fuel 
Cell stack 
Agilent Electronic Load 
= 


Fig. 2. Test bench. 
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be referred as “fresh” while the latter case will be referred as 
“aged”. 


4. Results 
4.1. Experimental results 


Fig. 3a and b, respectively, show the dynamic responses of 
each cell of the fresh and the aged stack to a steep change of the 
current density from 0 to 0.5 Acm™?. The graphs exhibit three 
main parts: 


(i) During the first second, a rapid decrease of the poten- 
tial from the OCV (about 1.06 V, not shown on Fig. 3) 
is observed, followed by an overshoot behavior, with a 
transient maximum in the cell voltage. Such behavior was 
already reported in the literature [11,13,22—24]. In addition, 
the amplitude of this overshoot is magnified by the ageing 
process: 24+4mV per cell for the aged stack instead of 
10+4mV per cell for the fresh stack. 

(ii) In the range between 1s and about 15s, a continuous 
decrease of the potential is observed. 

(iii) After 15s, an almost steady-state plateau is observed. 
Beyond 20s, the cell voltage gradually increases due to 
the internal temperature increase due to the Joule effect 
and the exothermic electrochemical reactions. In addition, 
the period requires to reach the new steady-state is higher 
for the cell located at the gases outlet than for the other 
cells: 20s instead of 15s. This time strongly depends on 


(a) 0,65 


© Cell #1 (gases outlet) 
o Cell #2 
0,5 


© Cell #1 (gases outlet) 
o Cell #2 


0,35 


15 


t(s) 


Fig. 3. Evolution of each cell voltage after a current step from OCV: (a) fresh 
stack and (b) aged stack. 
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Fig. 4. Evolution of each cell voltage during the first 60 ms after a current step 
from OCV: (a) fresh stack and (b) aged stack. 


the stack’s history and appears to be the most critical for 
the stack’s performances. 


It is worth mentioning that the cell at the gases inlet (cell 
#5) exhibits the best performances while the one at the gases 
outlet (cell #1) exhibits the lowest performances. So, the voltage 
difference between the first and the last cell into the stack is about 
30 mV for the fresh stack and about 80 mV for the aged stack. 

For further investigation, we magnified these evolutions in the 
first 60 ms after the current step on Fig. 4a and b. This graphs 
show that the cell voltage reaches a first plateau about 30—40 ms 
after current switch on. In addition, despite the short time cells 
performances are strongly decreased by the ageing, correspond- 
ing response time is similar for both the fresh and the aged 
stack. 

The data presented in Fig. 4a and b allows us to determine the 
cells ohmic resistance. Indeed, when an electrochemical cell is 
submitted to a current step from zero, its voltage first exhibits a 
steep change due to the ohmic drop. Then, it evolves as Faradaic 
processes and double-layer capacitance charge/discharge pro- 
ceed simultaneously, until the latter phenomenon ends [30]. As 
a consequence, if a fast acquisition system like the one used 
in this study monitors the voltage evolution, the voltage dif- 
ference between the last measurement before the current setup 
and the first measurement after the current setup can provide a 
good approximation of the cell’s ohmic resistance. The shorter 
the time duration between these measurements is, the better the 
approximation is. From the data presented in Fig. 4a and b, the 
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4 Experiment : cell nr 3 


aaa 


Fig. 5. Evolution of cell voltage during the first 60 ms after the current step, 
for a fresh stack versus model fitted for Cy =254F m7, Jeq = 110 A m~? and 
so =0.824. 


cell’s internal resistance can thus be calculated as follow: 


U — U= 
Ros ON as, (16) 
JimpSe 
Obtained values for fresh and aged cells were, respectively: 
1.05 +0.05 and 1.25+0.05mQ. These values were used as 


input parameters in our model. 


4.2. Simulation results 


Only three parameters were used to fit the evolutions of the 
simulated transient response to the experimental results: the ini- 
tial water volume fraction so, the equilibrium current density 
jeq and the double-layer capacitance Cy. Based on the previ- 
ous PEMFC dynamic model, a procedure is performed to fit the 
experimental results obtained for the middle cell (cell #3) of both 
fresh and aged stack. This cell was chosen because thermal and 
fluidic edge effects are minimized. The fitting was obtained by a 
Newton-Raphson algorithm that solves non-linear least square 
problems, using the Matlab optimization toolbox (LSQCURVE- 
FIT function). For the fitting, the data points corresponding to 
the overshoot behavior were removed, since our model does not 
take this phenomenon into account. For consistency, the exper- 
imental data points used for the curve fitting in the range 1—10 s 
were decreased by the magnitude of the observed overshoot. 
In addition, the experimental data points beyond 10s were not 
used for the curve fitting because our model does not take into 
account the cell temperature increase due to the heat production 
by the stack that impacts the cell voltage evolution. 

Figs. 5 and 6 compare experimental and simulated cell volt- 
age evolutions for the fresh stack, respectively, in the first 
60 ms and 10s after the current step. A very good agreement 
between experimental and simulated voltage evolutions were 
obtained for the following set of parameters: Cg) = 254Fm7?, 
Jeq = 110 A m7? and sọ =0.824. The observed slight discrepan- 
cies can be directly associated with the model assumptions: 
overshoot behavior and temperature evolution were not taken 
into account. Interestingly, the corresponding value of so is about 
80%, indicating that a high amount of liquid water is already 
present in the GDL before starting the cell [31]. 


4 Experiment:cell nr3 


— Model 


t(s) 


Fig. 6. Evolution of cell voltage during the first 10s after the current step, 
for a fresh stack versus model fitted for Cy =254F m7, Jeq = 110Am~ and 
so =0.824. 


Fig. 7. Evolution of cell voltage during the first 60 ms after the current step, 
2 


for a aged stack versus model fitted for Ca =135Fm™, jeq = 10Am~? and 
so =0.845. 


Figs. 7 and 8 compare experimental and simulated cell volt- 
age evolutions for the aged stack, respectively, in the first 60 ms 
and 10s after the current step. The model provides a good pre- 
diction of the cell’s time response and voltage level for the 
following set of parameters: Cg) =135F m7, Jeq = 10 A m~? 


o Experiment : Cell nr3 


Model 
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Fig. 8. Evolution of cell voltage during the first 10s after the current step, 
2 


for a aged stack versus model fitted for Ca =135Fm™, jeq = 10 A m~? and 
so =0.845. 


356 Ph. Mocotéguy et al. / Journal of Power Sources 167 (2007) 349-357 


and sg =0.845. A higher value of the initial liquid water con- 
tent so is required to fit experimental cell potential change in 
the case of the aged stack. This means that water management 
becomes more and more drastic to maintain acceptable perfor- 
mances. These results are in accordance with those obtained 
by Pasaogullari and Wang [32] who explained the decrease of 
PEMFC performances by the formation of liquid water in the 
cathode. This liquid water in the pore walls may reduce the effec- 
tive pore size available for gas motion leading to an increase of 
the mass transport limitations. 

The double-layer capacitance values (for the all electrode sur- 
face: 300 cm?) obtained from these fitting were 7.62 and 4.06 F 
for the fresh and the aged stack, respectively. Such decrease 
in the double-layer capacitance seems to indicate that the AL 
is affected by the ageing process. From these values, we can 
estimate the very minimum time necessary to charge the double- 
layer capacitance assuming that all the applied current is used 
for that. Indeed, as previously explained, after the steep change 
due to the ohmic drop, the voltage evolves as Faradaic pro- 
cesses and double-layer capacitance charge/discharge proceed 
simultaneously, until the latter phenomenon ends and the volt- 
age stabilizes. Thus, by setting js =0 in Eq. (13) and using 
experimental values from Figs. 5 and 7, the values of the very 
minimum time needed to charge the double-layer capacitance 
can be calculated as follow: 


dt = ——(Uo5 ms — U0.06s) (17) 


Obtained values were 15.3 and 10.9 ms for the fresh and the 
aged stack, respectively. This validates at posteriori the time 
limit of 0.1 s we have chosen in our calculation, beyond which 
the double-layer capacitance was considered as charged. 

Fig. 9 compares the evolutions of the calculated Faradaic and 
capacitive currents, for the aged and the fresh stack. This graph 
clearly shows that the double-layer capacitance is charged faster 
in the aged stack than in the fresh stack. But, in both cases, the 
time needed to charge this double-layer capacitance is lower than 
40 ms, and this is consistent with the time limit of 0.1 s chosen 
in our calculation. 
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Fig. 9. Comparison between the simulated evolutions of the Faradaic and the 
capacitive currents of the fresh and the aged stacks. 
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Fig. 10. Comparison between the simulated evolutions of the oxygen interfacial 
pressure in the fresh and the aged stacks. 


Fig. 10 compares the evolutions of the simulated oxygen pres- 
sure at the GDL/electrode interface for both fresh and aged stack. 
We notice that an induction period of about 10 ms is predicted by 
the model as the double-layer capacitance has begun to charge. 
Then, oxygen interfacial pressure decreases as an increasing 
part of the current is used to perform the electrochemical reac- 
tion. Then, the interfacial pressure evolution exhibits another 
time constant at about 150 ms, after which the interfacial partial 
pressure has decreased by 40% in the case of the fresh cell and 
by 51% in the case of the aged cell. This indicates that the oxy- 
gen reduction reaction rapidly becomes mass transfer limited. 
In addition, the model predicts that the electrode/GDL rapidly 
operates in starved oxygen conditions: the predicted interfa- 
cial pressure has decreased to 5 x 104 Pa after 7 and 19s of 
operation, for the aged and the fresh stack, respectively and to 
2.5 x 10* Pa after 14s of operation in the aged stack. In both 
cases, a linear decrease of the interfacial oxygen pressure with 
time is observed beyond 200 ms of operation. This corresponds 
to the linear increase of the liquid water volume fraction that can 
be calculated from Eq. (15) at a 0.165% s7! rate. 


5. Discussion 


The influence of the position of the cell in a stack operating 
in dead-end mode on its performances can be explained by two 
simultaneous phenomena: a water accumulation toward the out- 
let of the stack and a higher reactant gases partial pressure at 
the inlet of the stack. This is the consequence of the operating 
mode chosen for this study. Indeed, in dead-end mode, the gases 
are fed on one hand of the cell and their pressures are regulated 
on the basis of the value measured at stack’s inlet. And, what- 
ever the quality of the gas distribution system within the stack, 
there is a pressure gradient between the cells located at both 
ends of the stack. The fed gases are progressively consumed by 
the cells located at the stack’s upstream while produced water 
progressively accumulates at the stack’s downstream, pushed by 
the inlet gases velocity. This can explain why the cell located 
at gas inlet exhibits better performances than the one located at 
gas outlet. In addition, if we consider that the ageing process 
have caused a liquid water accumulation inside the stack, the 
decrease in each cell performances and the enhancement of the 
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performance differences can be easily explained. Indeed, as can 
be seen from Eqs. (3), (6) and (8), an increase in the water con- 
tent causes a decrease in the GDL porosity and thus a slower 
diffusion of the gases. This in turn increases the voltage losses 
by mass transfer limitation, as shown by our simulation results 
that indicate a strong and rapid decrease of the oxygen pressure 
at the GDL/electrode interface (Fig. 10). 

In addition, beyond a critical water content, the apparent 
exchange current density is also affected by water accumula- 
tion: the values of jeq that correctly fit our experimental results 
in the case of the aged stack is about 10 times lower than the 
one obtained by the fitting of the results obtained for the fresh 
stack. Such decrease in the actual exchange current density can 
be explained by a change in the accessibility of the active surface 
for the oxygen gas. Even if our simple model has not taken into 
account the complex structure of the active layer, it can be stated 
that liquid water partly fills the gas pore of the AL too. Diffusion 
path of dissolved oxygen in this layer is thus increased. In turn, 
diffusion limitation is increased, leading to a smallest effective 
equilibrium current density [31]. 

Finally, the very high values of the liquid water content 
(between 82 and 85%) that allow us to correctly fit our exper- 
imental results indicate that the water management issue in a 
fuel cell operating in dead-end mode is mainly associated with 
pore flooding while the risk of a membrane drying by insuffi- 
cient amount of water seems unlikely. As a consequence, the 
inlet gases humidification seems to be unnecessary to guaran- 
tee a good fuel cell performance during operation in dead-end 
mode. 


6. Conclusions 


The voltage evolution of a fuel cell operating in dead-end 
mode after a current step exhibits two different response times: 
one at about 40 ms, corresponding to the charge of the double- 
layer capacitance and the other one at about 15-20 s. The first 
time constant is independent of the stack’s history but the cor- 
responding voltage level is strongly impacted. An influence of 
the cell’s position in the stack on its performance has also been 
observed: the cell located at the gases inlet exhibits better per- 
formances than the one located at the dead-end. This difference 
was explained by a water accumulation at the outlet and by 
higher gases partial pressure at the cell located in gases inlet. 
The differences between the cell performances are enhanced by 
the ageing. 

The one-dimensional model we have developed is able to 
fit our experimental results with a good accuracy. The only 
observed discrepancies are linked with physical phenomena that 
were not taken into account in our model (increase in temper- 
ature by Joule effect and water distribution in the membrane) 
or to its simplification (the volumic structure of the active layer 
was not taken into account). Our model indicated that, even a 
fresh cell contains a high amount of liquid water when stored 
at room temperature and that the cell performance was rapidly 
limited by oxygen mass transport inside the GDL. Further exper- 
imental and simulation studies will be performed to improve the 


prediction’s accuracy of the model and to better understand the 
water management effect in a fuel cell operating in dead-end 
mode. 
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